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Abstract 

Sharkovskii proved that the existence of a periodic orbit in a one- 
dimensional dynamical system implies existence of infinitely many pe- 
riodic orbits. We obtain an analog of Sharkovskii's theorem for periodic 
orbits of shear homeomorphisms of the torus. This is done by obtain- 
ing a dynamical order relation on the set of simple orbits and simple 
pairs. We then use this order relation for a global analysis of a quantum 
chaotic physical system called the kicked accelerated particle. 

1. Introduction 

Given a dynamical system (X,f), a key question is which periodic 
orbits exist for this system. Since periodic orbits are in general difficult 
to compute, we would like to have the means to deduce their existence 
without having to actually compute them. 

Sharkovskii addressed the dynamics of continuous maps on the real 
line. He defined an order < on the natural numbers, Sharkovskii's or- 
der, and proved that the existence of a periodic orbit of a certain period 
p implies the existence of an orbit of any period q <p. We say the q 
orbit is forced by the p orbit. This offers the means of showing the 
existence of many orbits if one can find a single orbit of "large" period. 
For a dynamical system depending on a single parameter, if periodic 
orbits appear when we change the parameter, they must appear accord- 
ing to the Sharkovskii's order. Hence, Sharkovskii's theorem gives the 
global structure of the appearance of periodic orbits for one dimensional 
systems. Ever since the eighties there has been interest in obtaining 
analogs for Sharkovskii's theorem for two dimensional systems (see [I] 
and [IS]). 

A homeomorphism of a torus is said here to be of shear type if it 
is isotopic to one Dehn twist along a single closed curve. Let h be a 
shear homeomorphism, and let x be a periodic orbit of h. We can then 



define the rotation number of discussion in Section [2| Thus, a 

rational number in the unit interval [0, 1] is associated to each orbit. 

We consider orbits up to conjugation: orbits (x, f) and (y,g) are 
similar if there exists a diffeomorphism h of the torus T 2 which takes 
orbit x onto orbit y and hfh~ l is isotopic to g rel y. We define be- 
low a specific family of periodic orbits we call simple orbits. In this 
family there is a unique element up to similarity corresponding to each 
rotation number; hence they can be specified by their rotation num- 
bers. We emphasize it is not true in general that an orbit of a shear 
homeomorphism is characterized by its rotation number. 

These orbits are analyzed in Section [2} As it turns out, one simple 
orbit is indeed simple, and does not force enough dynamics. Hence 
we turn in Section [3] to analyze pairs of orbits. Two coexisting simple 
periodic orbits can form a simple pair and these are considered. The 
pairs do force some more interesting dynamics, as follows. For a pair 
of simple orbits of rotation numbers — and — to constitute a simple 
pair, it is necessary that the rotation numbers be Farey neighbors, i.e. 
|p2<Zi — Pi<?2| = 1- We denote such a pair of rational numbers by ^ V ^. 

We now define an order relation on the following set V of rational 
numbers and pairs in the unit interval, 

V = {q\q E Q D [0, 1]} U {q V p\q,p E Q D [0, 1]}. 

Define the order relation on V to be 

gVj)^r<^>r6 [q,p] 

and 

gVj)^rVs<^>r,sG [q, p) . 

Where we denote by [q,p] the interval between q and p, regardless of 
their order. 

Theorem 1.1. The order relation ^ on V describes the dynamical 
forcing of simple periodic orbits. Namely, the existence of a simple pair 
of periodic orbits with rotation numbers q\/p in V implies the existence 
of all simple orbits and simple orbit pairs of rotation numbers smaller 
than g V p according to this order relation . 

This is the main result of this paper and the proof is completed in 
Section |4j The idea of the proof is as follows. Consider the torus punc- 
tured on one or more periodic orbits of a homeomorphism h. Then h 
induces an action on this punctured torus, and on its (free) fundamen- 
tal group. Now apply the Bestvina-Handel algorithm to this dynamical 
system. The idea of using the Bestvina-Handel algorithm was used by 
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Boyland in [6] and he describes the general approach in [5]. In our 
case, after puncturing out a simple pair of orbits, applying the algo- 
rithm yields an isotopic homeomorphism which is pseudo-Anosov. The 
algorithm also offers a Markov partition for this system and we use 
the resulting symbolic representation to find that there are periodic 
orbits of each rotation number between the pair of numbers we started 
with. Then we directly analyze the structure of the pseudo-Anosov 
representative to show that all these orbits are in fact simple orbits. 
Furthermore any two of them corresponding to rotation numbers which 
are Farey neighbors form a simple pair. Finally we establish the iso- 
topy stability of these orbits using results of Asimov and Franks [2] and 
Hall [11 J - Thus, the orbits exist for any homeomorphism for which the 
simple pair exists, and are forced by it. 

One should compare this result with a very strong theorem of Doeff 
(see Theorem 3.6), where existence of two periodic orbits of different 
periods for a given shear homeomorphism h implies existence of peri- 
odic orbit of every intermediate rotation number. However an explicit 
description of these orbits is not given, while our theorem implies ex- 
istence of periodic orbits and pairs of orbits with precisely described 
behavior. 

This research was originally motivated by a question we were asked 
by Professor Shmuel Fishman: Is there a topological explanation for 
the structure of appearance of accelerator modes in the kicked particle 
system. In section [5] we give a description of the kicked particle system. 
This system turns out to be described precisely by a family of shear 
homeomorphisms of the torus. The existence of accelerator modes is 
equivalent to existence of periodic orbits. The global structure of this 



system is thus given entirely by the order relation in Theorem 1.1 while 



it cannot be directly computed due to the complexity of the system. 

The authors would like to thank Professor Shmuel Fishman for offer- 
ing valuable insights, Professor Italo Guannieri for some critical advice, 
and Professor Philip Boyland for many indispensable conversations. 



2. Simple orbits 

Let {xi, ,xn} be a set of points belonging to one or more pe- 
riodic orbits for a homeomorphism / of a surface S. The dynamical 
properties of this set of orbits are captured by the induced action of / 
on the complement S — S \ {xi, ....,xn} In a sense that will shortly 
become clear. Choose any graph G which is a deformation retract of 
the punctured surface S . A homeomorphism of S then induces a map 
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on G. The converse is also true: a given map of G determines a homeo- 
morphism of So up to isotopy. Therefore we specify the periodic orbits 
we analyze in terms of the action on a graph which is a deformation 
retract of the surface after puncturing out the orbit. 

Denote by Go the graph obtained by attaching N small loops to the 
standard unit circle at the points exp(^i), j = 0, . . . , N — 1. 

Definition 2.1. We call a periodic orbit x = {xi, . . . ,xn} for a shear 
homeomorphism / on the two-torus a simple orbit if the following hold. 
(1) There can be found a graph G which is a deformation retract of 
T = T 2 \ x as on Figure [l] such that G is homeomorphic toGg. 



Figure 1. A standard graph for a simple orbit 

We denote the loop in G corresponding to the unit circle the 
horizontal loop and the loops attached to it the vertical loops. 
(2) There exists a homeomorphism / of T isotopic to / rel x (i.e. 
the isotopy is fixed on x) so that a neighborhood of G is invari- 
ant under / and the induced action on G satisfies: (a) There 
exists a fixed number k e {0, .., iV — 1} such that each vertical 
loop is mapped k loops forwards (clockwise along the unit cir- 
cle) to another vertical loop, (b) The horizontal loop is mapped 
to itself with one twist around one of the vertical loops. 

Remark. A homeomorphism h for which we are given a simple periodic 
orbit must be of shear type as we can deduce from the action on the 
homology of the non-punctured torus. 

For a shear homeomorphism there exists a basis for the first homol- 
ogy for which the induced map is represented by the matrix 

From here on we refer to any two axes given by an homology basis that 
gives us the above representation as standard axes. The horizontal loop 
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FIGURE 2. The action on a standard graph for a simple orbit 



and one of the vertical loops in a graph for a simple orbits constitute 
a standard basis. 

Definition 2.2. Let h be a shear homeomorphism, and let h be a lift 
of h to the universal cover (a plane). For any periodic point x of h 
of period p, h p maps any lift x of x the same integer number q along 
the horizontal axis away from x, in a standard choice of axis (and x is 
possibly mapped some integer number along the vertical axis as well). 
We can then define the rotation number of x to be p(x) = | mod 1. 

The rotation number does not depend on the lift h of h. 

Remark. In the case of a homeomorphism isotopic to a Dehn twist on a 
torus, which is our interest here, it can be easily shown that the abelian 
Nielsen type equals exactly the rotation number defined above. 

There exists a simple orbit for any given rotation number r, and it 
is unique up to similarity. Denote the similarity class by f. In the 
following we use the word vertical to describe the y axis, in a standard 
choice of axis for / (the direction along which the twist is made). 

Lemma 2.3. Let x be a periodic orbit for a shear-type homeomorphism 
f of T 2 for which there exists a family of vertical loops such that they 
bound a set of annuli each containing one point of the periodic orbit, 
and this family is invariant under a homeomorphism f isotopic to f 
rel x. Then x is a simple orbit. 

Proof. Choose a vertical loop / of the invariant family. / is orientation 
preserving, and so is /. The first loop to the right of I is therefore 
mapped to the first loop to the right of f(l). Hence, the vertical loops 
in the invariant family are all mapped the same number of loops to the 
right. 

Now we have to find a horizontal line with the desired image. We 
write the invariant family of loops as {/.;}f =1 , where p = period(x), 



ordered along the horizontal axis. We choose another family of vertical 
loops {mj}f =1 , such that m 8 is contained in the annulus between U and 
h+i {lp+i — h), and passes through the periodic point Xi also contained 
in this annulus. Choose a point a\ / ii on mj. / can be adjusted in 
such a way that the new homeomorphism / leaves both families of 
vertical loops invariant, and in addition, so that a\ be a periodic point 
of / with period p. We denote the orbit of a% by {dj}f =1 where Oj G 
for 1 < % < p. 

Choose a line segment n\ connecting a\ to a 2 , so it crosses the annulus 
between m\ and m 2 from side to side. We choose rij + \ to be the line 
segment f^{n\) for 1 < j < p — 2. The boundary points of rij and n*, 
coincide whenever they lie on the same vertical loop. 

Now, we look at the horizontal loop n = (Jj=i n r Each segment of n 
is mapped exactly to the next segment, except n v ^\ which is mapped 
into the annulus between m\ and 1712. Since the mapping class group 
of an annulus is generated by a twist with respect to any loop going 
once around the annulus we may assume, that n p is mapped to rii plus 
a number of twists along such a loop. On the other hand we know that 
f(n) is homotopic to itself plus one twist in the negative direction (on 
the closed torus), so / maps n to itself plus one negative twist along 
this loop. By further adjustment of / we may assume the twist is made 
along I2. Thus the union of the vertical family with n chosen as 
above constitute a graph showing x to be a simple orbit. □ 

The Thurston-Nielsen classification theorem, see [7], states that any 
homeomorphism / on a compact connected oriented surface of negative 
Euler characteristic is isotopic to a homeomorphism / which is 

(1) pseudo Anosov, or 

(2) of finite order, or 

(3) reducible. 

where a homeomorphism <p is called of finite order if there exists a 
natural number n such that <p n = id. A homeomorphism <p is called 
pseudo-Anosov if there exists a real number A > and a pair of 
transverse measured foliations (J ru ,fi u ) and (JF s ,/i s ) with 0(jF u ,yu") = 
A/x u ) and 0(^-" s ,/i s ) = (J ru 1 \t^ u )- A homeomorphism on a sur- 
face M is called reducible if there exists a collection of pairwise disjoint 
simple closed curves V = {Ti, T^} in int(M) such that 4>(T) = V 
and each component of M \ T has a negative Euler characteristic. The 
representative / in the isotopy class of / which is of one of the three 
forms above is called the Thurston-Nielsen canonical form of /. 
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Of course homeomorphism / is reducible with respect to a simple 
orbit since it contains an invariant family of loops {/j} and the comple- 
ment of the invariant family consists of punctured annuli (which have 
negative Euler characteristic). 

Example 1. Not every periodic orbit for a shear homeomorphism is 
reducible. Consider the homeomorphism h described on Figure |3j It 
takes the graph on the left of Figure [3] to the graph on the right and is 
a shear homeomorphism. It has a periodic orbit of order 2, shown on 
the pictures, with rotation number 1/2 and it is pseudo-Anosov in the 
complement of the orbit. 




Figure 3. 



3. Simple orbit pairs 

Let x and y be two coexisting simple periodic orbits, for a homeo- 
morphism / of T 2 (/ must be of shear type), with rotation numbers ^ 
and H respectively. Assume p\ > P2, i.e., y has lesser period than x. 

Definition 3.1. We call the pair of orbits a simple pair if 

• We can find an embedded graph G in their complement home- 
omorphic to Go as on Figure |4| 

Each component in the complement of the graph is a topolog- 
ical rectangle which contains exactly one point of orbit x and 
at most one point of orbit y. 

• The homeomorphism / acts on this graph in the following way: 
each vertical loop except one moves to another vertical loop, 
there is one vertical loop denoted / such that f(l) is a vertical 
loop m plus a small loop around one point of the (shorter) y 
orbit, in a rectangle adjacent to line m on the right (as on Figure 
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Figure 4. 



[5]) or on the left, and the horizontal line is mapped to itself plus 
a twist in the negative direction around /(/), as on Figure [5] 




Figure 5. 



The graph which appears in Definition 3^ divides the torus T 2 into p\ 
rectangles. The homeomorphism / moves each vertical loop the same 
distance, say s rectangles, to the right except for the small additional 
loop for line I. Let Rq be the rectangle adjacent to m in which the 
small loop in the image of the graph occurs. Rq must contain exactly 
one point of each orbit. We denote these points x and y respectively. 
Under p\ iterations of / the point Xq runs q\ times around the whole 
torus, that is q\p\ rectangles to the right. So, pis = q±pi and s = q\. 

The point yo is mapped to itself after p 2 iterations. Under each iter- 
ation the image of yo is mapped s(= q\) rectangles to the right, except 
the last iteration under which it is moved an additional rectangle to the 
right or left. Altogether it has moved p 2 q± ± 1 rectangles. At the same 
time, it is mapped around the torus q 2 times, hence q 2 p± rectangles. 



This means £>2<Zi ± 1 = Thus for a simple pair of periodic orbits 

the rotation numbers — and — are Farey neighbors and the additional 

pi P2 JO 

loop is on the right (as on Figure IBT) if and only if ® > Denote 



P2 Pi 

by r V s~the similarity class of a simple pair corresponding to a pair of 
Farey neighbors r V s. 

Consider again the points xo and yo in the rectangle i?o- Continue 
the notation to all the points of x and y by Xi = T l (xo) and yi = T z (y ). 
We draw a small loop around each of the points of y. The union of 
these loops will be the peripheral subgraph P for the Bestvina-Handel 
algorithm, since we may assume the union of these loops to be f- 
invariant. Now we consider separately two cases. Case 1 will be the 
case in which m is the left boundary curve of Rq, while in case 2 it is 



21 
Pi 



< 



<J2_ 

P2 



and in the 



the right boundary (in other words in the first case 
second case ^ > ^). Choose some point on the loop around yo and 
connect it, by a curve l , to a point on the section of the horizontal line 
in R , in case 1 from below the segment and in case 2 from above. 




Figure 6. 



Then f(lo) is a curve connecting the loop around y± and the corre- 
sponding horizontal segment. We denote it by li, and do the same for 
each yi. After adding the above edges to the graph case 1 is topologi- 
cally as in Figure [6] 

The inclusion G So is a homotopy equivalence (where So is the 
punctured torus). We know the action of / on all edges of G except 
for the curve l qx -\ connecting y qi -i and the horizontal segment in the 
corresponding rectangle. It's image is a curve connecting the horizontal 
segment in the rectangle adjacent to m which is not Ro to the loop 
around y . This image might wind around a disk containing y and 
x qi -i as in Figure [7] 
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Figure 7. 



The graph and its image for case 2 are exactly the same except that 
the loops are connected to the horizontal segments from above. We 



shall prove in Proposition that we may assume that the image of 
the segment l qi -i has no winding. Hence we draw from now on the the 
graph images without winding, and we may assume the graphs given 
in Figure [8] also have an invariant neighborhood by a further isotopy 
of /. 
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Figure 8. The standard graph for a simple pair, case 
1 on the left and case 2 on the right. 



The action of / (up to isotopy) on this graph is given by one of the 
actions on Figure [9j drawn in some regular neighborhood of the graph, 
where each vertical loop moves qi loops to the right. 

Proposition 3.2. Let {x,y} be a simple pair with the graph as on 
Figure^and with windings as on Figure^ Then there may be chosen a 
different invariant graph, which also makes {x, y} a simple pair, whose 
image is without winding. 
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Figure 9. The action on a standard graph for a simple 
pair, for both cases respectively. 



Proof. To simplify the picture we prove the proposition for rotation 
numbers | and | . The general proof proceeds in the same way. We 
start by looking at a simple pair {x, y} of rotation numbers | and | for 
a homeomorphism / with the corresponding invariant graph G given 
so that the action on it is without any twists, as on the left side of Fig- 
ures [8] and |9j We now choose a different system of curves (a different 
graph), in a small neighborhood of G as in Figure 10, which will serve 
as a new graph for the pair. 




Figure 10. 



Solid lines are the new vertical loops and dashed lines are the new di- 
agonal segments like in figure [8} The horizontal loop consists of the 
dashed lines and the long pieces of the solid lines. To move from left 
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to right along the horizontal line, move along a dashed line and turn 
to the left when meeting a solid line. Continue up along a vertical loop 
and then along the next dashed line. We add to this graph the periph- 
eral subgraph and the connecting segments and get the graph H as on 



Figure 11 It is clear that topologically the graph H has the same form 



as the graph on Figure M and that it has an invariant neighborhood. 




Figure 11. 



The reader can check (using the precise knowledge of the image of 
each edge of the original graph) that the action of / on the graph H 
has the properties required from a simple pair. Each vertical loop is 
mapped onto another vertical loop except for one loop I for which f(l) 
is equal to a loop m plus a loop around the next periodic point y p 
of the shorter orbit. The horizontal loop is mapped onto itself plus a 
negative Dehn twist along /(/). Consider the image s of the segment 
which connects the horizontal loop to the periodic point y v -\. When 
the action has no twist then s moves along the horizontal loop in its 
positive direction until it meets the original segment connecting to y p 
and then it follows along the segment. However in our case s goes first 
backwards along the horizontal loop than moves in the counterclockwise 
direction along the boundary of the "rectangle" adjacent to the vertical 
loop m and finally follows the horizontal loop and the segment to y p . 
This means that the action / on the graph H has one positive twist. 

We proved that a simple pair for a shear homeomorphism with a 
given graph and a given action without twists can be given another 
graph which also describes it as a simple pair and the action on the new 
graph has one positive twist. This process is reversible. Therefore, by 

12 



induction, we can add or remove any number of twists using a suitable 
graph. This implies Proposition 4.2. □ 



Hence for a simple pair the action on a spine is given by Figures 
|8] and [9j We can now apply the Bestvina-Handel algorithm (see [3]), 
endowing a neighborhood of G with a fibered structure in the natural 
way. The algorithm specifies a finite number of steps which we apply 
to the graph G, altering G together with the induced action on it, but 
without changing the isotopy class of / on T 2 \ (x U y). When the 
algorithm terminates, it gives a new homeomorphism / which is the 
Thurston Nielsen canonical form of /. 

For simple pairs, the action in each of the two cases above is easily 
seen to be tight, as no edge backtracks and for every vertex there are 
two edges whose images emanate in different directions. The action 
has no invariant non-trivial forest or nontrivial invariant subgraph and 
the graphs have no valence 1 or 2 vertices. This is the definition in [3] 
for an irreducible map on a graph. 

Definition 3.3. Assuming g, the induced map on the graph itself, does 
not collapse any edges, there is an induced map Dg, the derivative of 
g, defined on 

C := J \{{v, e)\v is a vertex of G, e is an oriented edge emanating from 

by Dg(v,a) = (g(v),b) where b is the first edge in the edge path g(a) 
which emanates from g(v). 

Definition 3.4. We say two elements (v, a) and (v, b) in L corresponding 
to the same vertex v are equivalent if they are mapped to the same 
element under D(g n ) for some natural n. The equivalence classes are 
called gates 



The gates in each of the cases above are given by Figure 12 , indicated 
there by small arcs. There is no edge which g sends to an edge path 
which passes through one of the gates - enters the junction through one 
arm of the gate and exits through the other. Such an irreducible map 
is efficient, i.e., this is an end point of the algorithm. Now, since there 
are edges mapped to an edge path longer than one edge, we arrive at 
our next theorem. 

Theorem 3.5. A homeomorphism f of the two torus for which a simple 
pair of periodic orbits exists is isotopic to a pseudo-Anosov homeomor- 
phism relative to this pair of orbits. 
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Figure 12. 



Let / be a shear type homeomorphism of the torus, and fix a lift / 
of /. Define the lift rotation number of a point x G T 2 to be 



p(x, f) = lim n . 



(f n (x)-x), 



n 



for any lift x of x, when the limit exists, where the subscript 1 denotes 
the projection to the horizontal axis. Define the rotation set p(f) of / 
to be the set of accumulation points of 



X U 



x e R and iigN 



n 



Then, the above theorem follows from the following much more gen- 
eral theorem by Doeff, see [9] and |10j . 

Theorem 3.6. (Doeff) Let h be a shear type homeomorphism ofT 2 , 
and fix a lift h of h. If h has two periodic points x and y with p(x, h) ^ 
p(y, h) then h is pseudo-Anosov relative to x and y. Furthermore, the 
closure of the rotation set is a compact interval, and any rational point 
r in the interior of this interval corresponds to a periodic point x G T 2 
with p(x, h) = r. 

In particular, Doeff proves existence of two orbits of different rotation 
numbers implies existence of an orbit for any rotation number between 
these two. But he does not give any characterization of these orbits. 
Example [T] shows two different orbits, both with rotation number equal 
to 1/2, one of which is pseudo-Anosov, and the other reducible. Thus 
the rotation number does not give much information about the orbit 
and in this sense this theorem does not give a satisfactory dynamical 
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understanding of what is happening in regions of coexistence of orbits. 
In contrast with Doeff's general theorem, we get results for a very 
specific family of periodic orbits, but for this family we are able to give 
exactly the orbits forced by others, as we show in section |4j 

In our case the canonical form / of / we get by applying the Bestvina- 
Handel algorithm is a pseudo Anosov homeomorphism. When this is 
the case, the algorithm gives a canonical way of endowing a regular 
neighborhood S of G with a rectangle decomposition {Ri, R N }. 
The decomposition is a Markov partition for the homeomorphism /. 

A Markov partition for a dynamical system offers a symbolic repre- 
sentation for the system in the following way. Let T, N be the subset of 
the full iV-shift (the set of bi-infinite series on N symbols), where iV is 
the number of rectangles in the decomposition. Let E be a subset of 
defined by 

S = {s = (..., s n , s n+1 , ....): R Sn H f-'Rs^ ± 0} 

On S we naturally define a dynamical system with the operator of the 
right shift denoted by a, and (£, a) is called the subshift corresponding 
to the dynamical system. E can be completely described by stating 
which transitions k — > m for k, m G {1,....,N} are allowed (i.e., for 
which k, m, f~ l R m H Rk 7^ 0). See [lj for the definitions and for a 
proof that in this case we can define a map ix : S — > So by 

vr : Pin=of n Rs_ n n....nf-"R Sn 
which satisfies the following properties: 

• 7T<7 = fir, 

• 7r is continuous, 

• 71 is onto. 

We take here the set of sequences with the Tichonoff topology. Thus 
a periodic point in the symbolic dynamical system which is just a peri- 
odic sequence corresponds to a periodic point in the original dynamical 
system. 

To obtain the Markov partition in our case as in [3J, we thicken the 
edges of the graph to rectangles. In particular, the rectangles can be 
glued directly to each other without any junctions. This can be done in 
a smooth way, endowing So with a compact metric space structure by 
giving a length and width to each rectangle, consistently. Each edge of 
the standard graph for the pair (figure [8]) corresponds to one rectangle, 
except the edge which is mapped to the loop around %/q. This edge we 
divide in two (this is necessary to avoid having a rectangle intersecting 
twice an inverse image of another rectangle). Now we have edges of 7 
different types on the graph. The vertical loops of the graph consist of 
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long edges we denote as A edges, and short edges we call B edges. The 
loops around the points of the y orbit and vertical segments connect- 
ing the loops to the diagonal edges we call C's and -D's respectively. 
In rectangles which contain two punctures and therefore two diagonal 
edges we call the upper ones L edges and the lower ones K edges in 
the first case, and the lower ones L edges, upper ones K edges, in the 
second case. The last type of edges are diagonals of once punctured 
rectangles, these we call M edges. 
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Figure 13. 



Next, we label the rectangles in order to have explicitly the transition 
rules: 

• For < i < p 2 — 1 denote the rectangle corresponding to the 
D edge connecting the loop around yi to the diagonal by r i+ i. 
Denote the rectangle corresponding to the C edge which is the 
loop around yi by r P2+i+ i. 

• For 1 < i < p 2 , denote the rectangles corresponding to the L 
and K edges connected to r\ by r 2p2+ i and r 3p2+i respectively. 

• Denote the rectangle corresponding to the A edge belonging to 
the vertical line we referred to as m by r 4p2+1 and the B edge 
which is part of the same line m as r 4p2+Pl+ i. 

• For the vertical line f l (m) denote it's A and B rectangles by 
r Ap2+1+i and r 4p2+Pl+1+i respectively for all 1 < % < pi - 2 

• For the vertical line / P1_1 (m), denote it's A edge as by r^ P2+Pl . 
There are two rectangles corresponding to the B edge as ex- 
plained above, denote the lower one by r Ap2+2pi and the upper 
one by r 4p2+2pi+1 . 
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• Label the pi — P2 remaining rectangles corresponding to the M 
edges by starting with the first of these to the right of m, and 
then continuing by the order along the horizontal axis, denoting 
them by r 4p2+ 2p 1+ 2, ••• ,r 3pi+ 3p 2+ i 

Finally, we can look at the diagram in figure 14, showing the set of 
rectangles and transitions in this Markov partition which we now use. 




Figure 14. Some of the rectangles in the Markov par- 
tition, where the arrows denote allowed transitions be- 
tween them 

A periodic symbolic sequence of allowed transitions gives as ex- 
plained a periodic point in the original dynamical system. Therefore 
by this diagram we can easily find other periodic orbits on the torus 
that must exist for /. We will later prove that these orbits are in 
fact simple, but this will require some more work. Hence, by this 
diagram we prove only existence of orbits with specified rotation num- 
bers. For every pair (n, m) of natural numbers, n, m ^ 0, by starting 
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from r pi+ 4 P2+ i, going n times around the first loop in the diagram 
{r Pl+4p2+ i, . . . , r 2pi+ 4j, 2 }, then going m-1 times around the second loop 
{r±, . . . ,r P2 } (and skipping it if m=l) and then returning through the 
final sequence {r2 P2 +i, . . . , r 3p2 } to r Pl+ip2+ i, we get a periodic symbolic 
allowed sequence, and so a new periodic orbit we denote O n ,m- These 
symbolic sequences are all different and hence so are the periodic orbits. 
We look at a point p G O n>m such that p is in the rectangle r pi+4p2+1 . 
For the first n ■ p\ iterations of p, corresponding to each time the upper 
loop in the diagram appears in the symbolic sequence of p, the images 
are contained in the B edges. The vertical loops are mapped under / 
retaining the same "horizontal distance" from the periodic points from 
the x orbit to their left. So, p is mapped a total distance of n • qi along 
the horizontal axis under f' 1 *. 

Similarly, point f npi (p), which lies in the rectangle T\ corresponding 
to D edge, is mapped a distance g 2 along the horizontal axis under each 
iteration of f P2 , for every occurrence of the second loop in the symbolic 
sequence of p. This is because the D edges retain their distance from 
the y orbit points below them. The final sequence in the symbolic 
representation of p until the return to the first loop also corresponds 
to the horizontal distance q±. These last points of the periodic orbit lie 
in rectangles corresponding to L edges. So p is mapped a horizontal 
distance of nqi + mg 2 under f n P± +m P2 . Hence, the new orbit O n ,m has 
rotation number " gl "|" mg2 . 

npi+mp2 

See [12] for a proof that any two Farey neighbors span this way all 
rational numbers between them, that is all rationals between — and — 

' Pi P2 

are of the form ragl "|" mi?2 . So we found a periodic orbit of any rational 
number between the two original rotation numbers — and — . 

We have found these periodic orbit to the canonical form of the 
homeomorphism we started with. We must now relate these to the 
periodic orbits of the original homeomorphism. Recall the following 
definition from [2]. 

A periodic point xq G S of period p for homeomorphism /q is called 
unremovable if for each given homomorphism fx with f t : fo ~ j\ there 
is a periodic point X\ of period p for /^and an arc 7 : [0, 1] — > S with 
7(0) = x , 7(1) = x\ and j(t) is a periodic point of period p for f t . 

It was proven by Asimov and Franks in [2j that every periodic orbit of 
a pseudo-Anosov diffeomorphism is unremovable. Thus orbits found for 
the pseudo-Anosov representative exist for any other homeomorphism 
in its isotopy class. This yields all these periodic orbits exist for the 
original homeomorphism / as well. Thus we get theorem 3^ for our 
specific case: 
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Theorem 3.7. If there exists a simple pair of orbits for a homeomor- 
phism f of the torus of abelian Nielsen types q and p which are Farey 
neighbors, there exists a periodic orbit for f with abelian Nielsen type 
equal to r for every rational number r between q and p. 



4. The order relation 



For any simple pair 



Pi 



V 
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l->2 ' 



The orbit On out of the family of 



new orbits we constructed above has rotation number equal exactly to 
91+92 This orbit corresponds to the symbolic sequence r f 



P1+P2 
r Pi+4p 2 +2 

diagram in figure [14 



'pi+4p 2 +l 



?"2pi+4j>2 — * r 2p 2 +l 



T 3P2 



So we have a list of rectangles, each containing 
exactly one periodic point from the new orbit 0\ \. We denote the point 
of Oi l that is in a rectangle rj by Oj. Graphically, assuming the first 
case map, when we draw the rectangle decomposition corresponding to 



the standard graph as in figure p\ we get Figure 15 




Figure 15. The intermediate orbit denoted by black 
circles. The gray areas are junction, which can be 
deleted, and the rectangles can be glued directly to one 
another. 



We will now show that for any simple pair ^ V the orbit Oi.i 
is a simple orbit, and forms a simple pair with each periodic orbit of 
the pair, that is with — and — . For the first assertion, we define 

Pl P2 ' 

a family of vertical loops as follows: we choose a vertical loop that 
crosses both rectangles corresponding to the m line and passes to the 
right of the periodic point o pl+ip2+ i. Denote this loop by A. It is shown 
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graphically in figure 16 All its images under / until the pist iteration 
are exactly of the same form, as the rectangles are simply mapped to 
the right without changing their forms. Its pist image is the first time 
it returns to the same rectangles, and is determined by the images of 
the corresponding vertical edges of the graph. These images are shown 
in figure [9j We use the fact / preserves orientation to determine the 
relation between the image of the curve and the points of On. We 



denote this image by B. It is shown in figure 16 




Figure 16. The intermediate orbit with the family of 
vertical loops 

By similar considerations, knowing the rectangles containing B in the 
original picture (Figures [8] and [9]) the rectangle adjacent to f pi ~ l {m) 
on the left contains a point of the y-orbit and therefore contains a 
rectangle of type L of Markov partition. This rectangle contains the 
point 03 P2 of the new orbit. Line A lies to the right of m and of o Pl+ 4 P2+ i 
therefore f Pl ~ l (A) lies to the right of o 3p2 and to the right of 02 P i+4 P2 - 
It follows that the line B = f Pl (A) lies to the right of o pl+4p2+ i and to 
the right of o 2p2 +i, as shown on Figure 16 Also B can be isotoped to 
the right of A relative to the points of the new periodic orbit. Next 
P2 — 1 iterations of / translate A and B and whole rectangle adjacent 
to m on the right to the right. We arrive at the rectangle adjacent to 
I on the left containing point o^ P2 of the new orbit. The point c>2 Pl +4 P2 
lies in a rectangle on the line /. The loop f P2 ~ 1 {B) lies to the right of 
o 3p2 and to the left of o 2pi+ 4 P2 therefore the loop C = f P2 (B) lies to the 



right of o pi+ 4 P2+ i and to the left of o 2p2 +i, as shown on Figure 16 Also 
p 2 iterations of / take line m to a distance p 2 <?i = P1Q2 — 1 rectangles 
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to the right, which means one rectangle to the left of line m. Since A 
and B are to the right of the point o Pl+ 4 P2+ i in line m and since this 
point moves to the leftmost point in the new periodic orbit shown on 



Figure 16, loop C must be to the right of it, as in Figure 16 The point 
c»2pi+i may be above or below the loop C but this does not change the 
discussion bellow. 

Note that if we disregard the orbits x and y of the original pair we 
can isotop C to A relative the points of the new orbit. This shows that 
the new orbit is a simple periodic orbit of length pi + p 2 , by Lemma 



2.3. Now we fill in the x orbit (the longer orbit) and consider a torus 
punctured at the y orbit and the new orbit together. We have the 
family of vertical loops f l (A) and the action on it is exactly as in the 
condition for a simple pair as the loop C can be isotoped to A plus a 
loop around y$. We choose a horizontal loop as the loop D on Figure 
17 Then its image D' is as shown on Figure 17 The image has the 
required properties. The orbit y together with the new orbit form a 
simple pair for the homeomorphism /. 




Figure 17. 

Next we fill in the y orbit and leave punctures at the x orbit and 
the new orbit. We choose the initial vertical loop A differently, as 
in Figure 18 This loop A is one rectangle to the right of m plus a 
loop on the left. After k = p 2 — 1 iterations of / it will move to 
line / which is q\ rectangles to the left of m. Indeed it will move to 
Qi(P2 ~ 1) + 1 — P1Q2 — Qi rectangles to the right of m which means qi 
rectangles to the left. The loop f k (A) looks like the loop A and lies 
to the left of the point r 3p2 and to the left of the point o 2pi+ 4 P2 . Next 
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iteration of / takes it to a curve which looks like /(/) but lies to the 
left of o Pl+ 4 P2+ i and to the left of 02 P2 +i- Since we filled the point yo we 
can isotop this loop to a vertical loop near m, which passes to the left 
of o pl+ 4P 2+ \. Subsequent iterations translate it to the right and f pl (A) 



is equal to curve B on Figure 18 




Figure 18. 

Next P2 — 1 iterations will take B to the loop near / which lies to the 
left of o 2 p 1+ 4p 2 . Next iteration of / takes this loop to a loop similar to 
/(£), but lies to the left of o 2p2+ i. Since the points of ?/-orbit are filled 



we can isotop it to the loop C on Figure [18] It can be further isotoped, 
relative to the x-orbit and the new orbit, to the loop A plus a small 
loop around xq to the left of A. If we choose the same horizontal loop 
as in the previous case , with the same image as before, we get the 
required action of / for a simple pair consisting of the x-orbit and the 
new orbit. 

Now we can continue by the same analysis for each of these two 
simple pairs, finding their Farey intermediate to be a simple orbit as 
well that forms a simple pair with each of them, and so on. It remains 
to prove the persistence of all these simple pairs under isotopies. For 
this, recall The following theorem from [TT] . 

Theorem 4.1. (Hall) Let S be a closed surface and let A be a finite 
subset of S. Let f be a homeomorphism of S which leaves A invari- 
ant. Let p = x\, . . . ,Xk be a finite collection of periodic points for f 
which are essential, uncollapsible, mutually non- equivalent and non- 
equivalent to points of A. Then the collection p is unremovable, which 
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means that for every homeomorphism g isotopic to f rel A there exists 
an isotopy f t rel A and paths Xi(t) in S such that fo = f , fi = g, 
Xj(0) = Xi, Xi(t) is a periodic point of ft of period equal exactly to the 
period of Xi. 

(This theorem is a generalization of the main result of Asimov and 
Franks in [2] to several periodic orbits. In fact this generalization was 
mentioned in [2] as a remark with a hint of a proof.) 

Recall also that if / is pseudo-Anosov in the complement of A then 
it is condensed and by [6] Lemma 1 and Theorem 2.4 each periodic 
point is uncollapsible and essential and points from different orbits are 
non-equivalent and points disjoint from A are not equivalent to points 
of A. 

Corollary 4.2. : Let T be a torus and let A be a finite subset of T. 
Let f be a shear-type homeomorphism ofT which is pseudo-Anosov in 
the complement of A. Let g be a homeomorphism ofT isotopic to f in 
the complement of A. If x is a simple periodic orbit for f then there 
exists a simple periodic orbit z for g with p(z) = p(x). If x,y is a 
simple pair of periodic orbits for f , one or both disjoint from A, then 
there exists a simple pair of periodic orbits z,w for g with p(z) = p(x) 
and p(w) = p(y). 

Proof. Chose points x\ and y\ from the orbits x and y . By Theorem 



4.1 there exists an isotopy f t and paths Xi(t) and yi(t) such that 2i(0) = 
x ii 2/i (0) = y%i Xiit) is a periodic point of f t of a fixed order p for all 
t and yi(t) is a periodic point of f t of a fixed order q for all t and 
2/i (t) = 2/1 (0) for all t if y(0) G A. For a given t all points in the orbits 
of xi(t) and yi(t) for / are distinct, they form a braid with p+q strands. 
They move when t changes and their movement can be extended to an 
ambient isotopy h t which is fixed on A. Then h t (f l (xi(0)) = f{(xi(0)) 
and h t (f l (yi(0)) = fl(yi(0))- Consider isotopy F t = h^ 1 f t h t . We have 
F t (f( Xl (0)) = f +1 (xi(0)) and F t (f( yi (0)) = f +1 (yi(0)) so F t is fixed 
on the orbits x and y. In particular x and y form a simple pair of 
periodic orbits for F\ (or x forms a simple periodic orbit for F\ if there 
is no y). But F\ = K{ x gh\ so h\{x) and h\(y) form a simple pair of 
periodic orbits for g. 



This concludes the proof of Theorem 1.1 



5. Global analysis of the kicked accelerated particle 

SYSTEM 



The physical system called the kicked accelerated particle consists 
of particles that do not interact with one another. They are subject 
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to gravitation and so fall downwards, and are kicked by an electro- 
magnetic field, i.e., the electro magnetic field is turned on for a very 
short time once in a fixed time interval. This electromagnetic field is a 
sine function of the height of the particle, hence the particles are kicked 
upwards or downwards by different amounts, depending on their po- 
sition at the time of a kick. For a short review of the results for this 
system see [5]. Experiments of this system were conducted by the Ox- 
ford group, see [TH], and the system was found to show a phenomena 
that is now called "quantum accelerator modes": as opposed to the 
natural expectation that particles fall with more or less the gravita- 
tional acceleration, it was found that a finite fraction of the particles 
fall with constant nonzero acceleration relative to gravity, as can be 



seen in Figure 19 




-200 -150 -100 -50 50 100 150 
Momentum 



Figure 19. Accelerator modes 
Experimental Data (taken from Oberthaler, Godun, d'Arcy, Summy 
and Burnett, see [16]) showing the number of atoms with specified 
momentum relative to the free falling frame as the system develops in 

time (the numbers on the y axis represents time by the number of 
kicks, while the z coordinate is proportional to the number of atoms) 



This is a truly quantum phenomenon having no counterpart in the 
classical dynamics. A theoretical explanation for this phenomenon was 
given by Fishman, Guanieri and Rebuzzini in [13] , and it establishes a 
correspondence between accelerator modes of the physical system, and 
periodic orbits of the classical map 
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(1) /:(^)-( J + ^+- /)+n ) modto 

Where the J coordinate corresponds to the particles momentum, and 
9 to its coordinate. This map is of shear type, and the acceleration for 
a periodic orbit with rotation number - is given by 



(2) a= 2 -^-Q 

V 

Hence, by analyzing the structure of existence of periodic orbits for 
the classical map above, we would be able to find which modes should 
be expected for which values of the parameters k and fl. We remark 
that actual experimental observation also requires stability of the pe- 
riodic orbits. It is important to stress here that since these parameters 
correspond to the kick strength and the time interval between kicks 
they can be controlled in the experiments as we wish, so results ob- 
tained for this system can be tested experimentally. When one plots 
the numerical results describing which periods exist for different values 



of k and fl one gets an extremely complicated figure, see figure 20 

An exact mathematical analysis of this system is extremely compli- 
cated. Perturbative methods have been used in [T3] to analyze the ex- 
istence of these "tongues" of periodic orbits in the region where k — » 0, 
as well as giving estimates on their widths. 

Look at the map / given by (|T| in regions where fl is equal 
for some rational number - in the unit interval, and small k. For a 
small enough k it can be seen both from the numerical results shown 



graphically in figure 20 and from perturbative arguments that in the 
above region a periodic orbit with period p exists. 

For small k the periodic points of this orbit must be pretty much 
equally spaced along the J axis, and we can choose (for k small enough) 
a family of vertical loops that are equally spaced at distance exactly fl 
apart, and each is at distance at least, say, 3k from any of the periodic 
points. 

The image of a loop parameterized by Ti(9) = ( ^° ] is given by 



J + ksin(J + 0) + fl \ / J + ksin(Q') + fl \ 

J + 6 ) = { 6' ) ^d so is very 

close (for small k) to another loop of the chosen family. It follows that 
there exists a map / isotopic to / rel the orbit which keeps this family 
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of curves invariant, and so, by lemma 2.3 all the periodic orbits seen 



in the tips of the tongues in Figure 20 are simple orbits. 



Note the rotation number of each of these orbits is equal exactly to 
the value of Q in the tip of the tongue (k = 0) as for very small k 
the J coordinate increases by an almost fixed value, close as we wish 
to Vt. And, by equation (|2| the rotation number | is related to the 
acceleration of the corresponding acceleration mode by 

a = 2 2vr - Q 
v 

So the topological meaningful numbers here are in fact also the ones 
with physical significance. While Q changes through the region in which 
this periodic orbit exists, | is of course a topological invariant and 
therefore fixed. Hence the acceleration vanishes on the line with fixed 
Q in the middle of each tongue, and changes signs when one crosses 
this line. This was measured experimentally in |14j . 

For any other point higher in the tongue which we can reach by an 
isotopy along which the periodic orbit exists, we also have the orbit 
is a simple orbit. We will assume, as is very natural and was checked 
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numerically for many cases, that the orbits remain simple throughout 
the region of each tongue. 

In some of the cases for which we drew a portrait of the phase space, 
we found that the fact the homeomorphism is isotopic to one which is 
reducible rel the periodic orbit is realized by the physical map itself, as 



seen in Figure 21 




Figure 21. Phase portrait for a two-orbit 
Drawn for k = j^2n and Q = n, the two-orbit which is clearly seen is 
a stable orbit with two stable neighborhoods drawn. There is another 
two-orbit present, at which the arrows point, and it is the stable and 
unstable manifolds for this unstable orbit which divide the phase 
space into non intersecting regions which do not mix. 



Here the phase space is truly divided into pieces. Each of the annuli 
in this decomposition is mapped to another, and returns to itself with 
one twist after p iterations of /. Therefore every periodic orbit must 
have a period which is a multiple of p. On the other hand, when an 
annulus is mapped to itself with one twist under an area preserving 
map (here under f p ), every rotation number in the unit interval exists 
for it (here we mean the standard annulus rotation number measuring 
the rotations around the annulus), and so every period exists, as for 
every rational number — there is a periodic point of order m which 
rotates n times around the annulus before it returns to itself. This 
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yields that for such a point in the parameter space, exactly all periods 
that are multiples of p exist. It is our belief that this situation is typical 
for the center of each tongue, that is for Q = 2tc^. At other points, 
namely in all point we have numerically checked outside the center of 
the tongue, orbits of coprime lengths may exist simultaneously. We 
believe that the coexisting orbits whose rotation numbers are Farey 



neighbors form a simple pair together, as in the example on Figure 22 



which shows a simple pair of orbits with rotation numbers 1/3 and 1/2 
found in the physical system. 




Figure 22. A pair of coexisting orbits in the physical system 
Drawn with a collection of curves on the torus and their images, 
which show this is a simple pair. 



This coexistence happens at a point in Figure 20 for which two 
tongues intersect. We assume that the same orbit persists through- 
out the tongue, and therefore we have at such a point two coexisting 
simple orbits. We believe that in all points of intersecting tongues 
coming from k = and Qi = Q2 = ^ which are Farey neighbors, 
Pi > P2, the coexisting orbits form a simple pair. 



Theorem 3/7 therefore implies that there are infinitely many periodic 
orbits for the parameters at a region of intersection of two such tongues, 
with rotation numbers equal to all rational numbers between the ones 
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of these two tongues. If we assume all these simple orbits present also 
come from tongues, this yields that each rational tongue between ^j- 
and H intersects each of these two tongues lower (along the k axis) than 
they intersect each other. In other words, following a path from a tip of 
a tongue upwards in the tongue, if it intersects a Farey neighbor tongue 
we know it intersects earlier all tongues of rational numbers between 



them. This determines the global structure appearing in Figure 20 of 
all accelerator modes in the physical system, as Sharkovskii's theorem 
determines it for one dimensional systems. 
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